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This  paper  is  concerned  with  the  stability  and  Robust  stabilization  prob¬ 
lem  for  2-D  continuous  systems  in  Roesser  model,  based  on  Generalized 
Kalman— Yakubovich— Popov  lemma  in  combination  with  frequency-partitioning  ap¬ 
proach.  Sufficient  conditions  of  stability  of  the  systems  are  formulated  via  linear  ma¬ 
trix  inequality  technique.  Finally,  numerical  examples  are  given  to  illustrate  the  effec¬ 
tiveness  of  the  proposed  method. 
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1.  INTRODUCTION 

Stability  of  2-D  continuous  systems  is  the  major  aim  in  all  researches,  in  order  to  guarantee  the  nor¬ 
mal  operation  of  systems.  In  relation  with  these  researches,  there  are  various  results  in  the  past  decades.  For 
example,  the  stability  of  2-D  continuous  systems  has  been  solved  lately  in  [1],  the  stability  margin  of  2-D  con¬ 
tinuous  systems  have  been  computed  with  a  new  method  in  [2],  FMI  based  stability  analysis  for  2-D  continuous 
systems  was  considered  in  [3],  Robust  stability  analysis  for  2-D  continuous-time  systems  was  obtained  in  [4], 
the  stability  analysis  based  on  the  quadratic  Fyapunov  function  was  obtained  in  [5],  iJoo  filtering  of  uncertain 
2-D  continuous  systems  with  time-varying  delays  was  considered  in  [6] .  In  addition,  the  Robust  stabilization 
and  control  design  have  been  studied  in  some  papers  as  well,  to  list  some  of  these,  authors  in  [7]  proposed  the 
robust  state  feedback  Hoo  control  for  uncertain  2-D  continuous  state  delayed  Roesser  systems.  Hoo  control 
of  2-D  continuous  switched  systems  have  been  investigated  in  [8],  multiobjective  H2/H00  control  design  was 
considered  in  [9],  FMI  based  robust  PID  control  has  been  solved  in  [10],  and  stabilization  of  two-dimensional 
continuous  systems  have  been  investigated  in  [11]. 

Recently,  attention  has  been  devoted  towards  the  Kalman— Yakubovich— Popov  (KYP)  lemma  in  [12], 
this  lemma  makes  equivalence  between  frequency  domain  inequality  (FDI)  characterizing  a  class  of  properties 
of  a  transfer  function,  and  a  linear  matrix  inequality  (FMI)  in  [13],  for  its  state  space  realization.  Therefore, 
authors  in  [14]  has  proposed  an  extension  of  the  KYP  lemma,  which  is  known  as  Generalized  KYP  (GKYP) 
lemma  for  the  case  of  finite  frequency  domain.  The  2-D  GKYP  lemma  is  obtained  for  Roesser  model  of  2-D 
continuous  systems  in  [15],  and  for  2-D  discrete  systems,  for  both  cases:  Fornasini-Marchesini  (FM)  and  for 
Roesser  models,  in  [16]  and  [17],  respectively.  The  GKYP  combined  with  the  frequency-partitioning  approach 
to  stability  analysis,  were  obtained  in  [18]  for  2-D  discrete  system,  and  for  hybrid  systems  in  [19,  20]. 

Motivated  by  the  Previous  research,  in  this  paper,  we  suggest  a  sufficient  conditions  of  stability  of 
2-D  continuous  Roesser  systems,  via  GKYP  lemma  and  frequency-partitioning  approach,  in  order  to  reduce 
the  conservativeness  of  the  existing  simple  2-D  continuous  Fyapunov  inequality.  Generally,  in  order  to  realize 
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a  series  of  novel  stability  conditions  for  our  system,  the  GKYP  lemma  is  applied  on  each  one  of  the  N  intervals 
of  the  entire  frequency  domain.  Moreover,  robust  stabilization  is  also  considered  based  on  the  proposed  stabil¬ 
ity  conditions.  Finally,  numerical  examples  are  given  to  demonstrate  the  effectiveness  of  the  proposed  method. 

Notation:  we  use  the  following  notation  throughout  this  paper.  The  superscript  T,  *,  -1  stand  for  matrix 
transpose,  matrix  complex  conjugate  transpose,  and  matrix  inverse,  respectively.  /  denotes  an  identity  matrix 
with  appropriate  dimension.  The  notation  P  >  0  (P  <  0)  means  that  matrix  P  is  positive  (negative)  definite. 
diag{.}  stands  for  the  block  diagonal  matrix,  ReX{.)  is  the  real  of  eigenvalue  of  a  square  matrix.  Matrices,  if 
their  dimensions  are  not  explicitly  stated,  are  assumed  to  be  compatible  for  algebraic  operations. 


2.  PROBLEM  FORMULATION  AND  PRELIMINARIES 

Consider  the  following  Roesser  model  for  2-D  continuous  systems  : 


dx^(ti,t2) 

dti 

9*2 


Ai  A2 
A3  A4 


X''{ti,t2) 


+ 


Bi 

B2 


U{ti,t2) 


(1) 


where  x^{ti,t2)  G  ^2)  G  K”’'  and  u{ti,t2)  G  are  the  horizontal  state,  vertical  state  and  input 

of  system,  respectively,  and  Ai,  A2,  A3,  A4,  Bi  and  B2,  are  real  matrices  with  appropriate  dimensions. 

For  real  numbers  ti  and  t2,  we  introduce  notations 


=  sup  11  X^{0,t2)  II,  x'^  =  sup  II  ||  . 

Assumption  1 

lim  II  a:(fi,0)  ||=  0  and  lim  ||  x{0,t2)  ||=  0. 

ti - yoo  t2 - ^00 

They  are  inferred  to  the  initial  condition  for  the  system  (1). 


In  the  stability  analysis  of  2-D  continuous  system  (1),  it  is  required  to  consider  the  zeros  of  the  2-D  characteristic 
polynomial  given  by 


<^(51,  S2)  =  det 


Sllnh  —  Ai  —A2 

—  A3  S2lnv  —  A4 


(2) 


It  is  known  in  the  literature  that  the  2-D  continuous  system  is  asymptotically  stable  if  and  only  if  67(81 ,  S2)  f  0 
V(si,  S2):  Pe(si)  >  0  and  Re{s2)  >  0. 

In  general,  this  condition  is  difficult  to  use  in  practice  to  verify  the  stability,  therefore,  another  method  will  be 
used  via  LMI. 


Lemma  1  Simple  necessary  conditions  for  asymptotic  stability  of  the  2-D  continuous  system  (1)  are  as 
follows: 

i)  Ai  is  Hurwitz  (i.e.  all  its  eigenvalues  have  negative  real  parts,  ReXi{Ai)  <  0;  i  =  1, ..,  n/i). 

ii)  A4  is  Hurwitz. 


Proof  From  (1)  for  A2  =  A3  =  A4  =  0,  we  obtain  the  state  equation  of  the  continuous  system  (  for  the 
fixed  0  <  ^2  G  P ) 


dx'^(ti,t2) 

dti 


=  AiX^{ti,t2) 


The  system  (3)  is  asymptotically  stable  if  the  matrix  Ai  is  Hurwitz. 
Similarly,  we  can  proof  ii). 


(3) 

□ 


Therefore,  we  assume  the  following  throughout  the  paper. 

Assumption  2  The  matrices  Ai  and  A4  are  Hurwitz. 

Lemma  2  Let  the  assumption  2  be  satisfied,  the  2-D  continuous  system  (1)  is  asymptotically  stable  if  and 
only  if 

S'(s)  =  ^3(3/ —  Ai)  ^A2 -f  A4,  s  =  jui  (4) 

is  Hurwitz  matrix  for  a;  G  M. 
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Proof  Letu{ti,t2) 

initial  condition,  we  get 


0,  and  taking  the  Laplace  transformation  of  system  (1)  for  ti  only,  and  under  zero 


sX^{s,t2) 

dX'’(s,t2) 

dt2 


Ai  A2 

■  X>^{s,t2)  ■ 

A3  A4 

.  X-{s,t2)  _ 

(5) 


Solving  (5),  we  obtain 

- ^  —  ^l)  ^  ^2  +  A4]X‘"  {s,t2) 


(6) 


System  (6)  can  be  regarded  as  a  1-D  continuous  system  with  complex  variable  s,  and  we  notice  that  the  variable 
t2  of  the  system  doesn’t  depend  on  the  variable  s.  So  the  1-D  continuous  system  (6)  is  asymptotically  stable  if 
and  only  if  [^3(5/  —  Ai)~^A2  +  A4]  is  Hurwitz  matrix  for  Re{s)  =  0.  Hence,  the  2-D  continuous  system  (1) 
is  asymptotically  stable  if  and  only  if  [^3(^0;/  —  Ai)~^  A2  +  A4]  is  Hurwitz  matrix  Vw  G  M.  □ 


Remark  1  Notice  that  when  interchanging  Ai  with  A4,  and  A2  with  A3,  the  2-D  continuous  system  (1)  is 
asymptotically  stable  if  A2  {jcol  —  ^4)“^  A3  -f  Ai  is  Hurwitz  matrix  \/uj  G  M. 


We  will  use  the  following  lemmas,  known  as  the  KYP  lemma,  the  GKYP,  and  the  Projection  Lemma,  respec¬ 
tively. 


Lemma  3  [12]  Let  matrices  A,  B,  and  0  =  0^  be  given,  if  det{jwl — A)  f  0  Vw  G  M.  Then  the  following 
two  statements  are  equivalent. 

( i)  For  any  w  G  K  U  00, 


{jwI-A)-^B 

I 


{jwl  -A)-^B 
I 


<  0 


(ii)  There  exists  a  symmetric  matrix  P  such  that 


A  B 
I  0 


0  P 
P  0 


A  B 
I  0 


-f  0  <  0 


(V) 


(8) 


Lemma  4  [14]  Let  the  matrices  0,  F,  $  and  4*  be  given,  and  denote  is  the  null  space  ofTi^jF,  where 
T(^  =  [  /  —joji  ].  The  inequality 


W*0iV^  <  0,  with  oj  G  [ijJi,u}2], 
holds  if  and  only  if  there  exist  Q  >  0  and  a  symmetric  matrix  P,  such  that 


(g)  P  -f  4*  (g)  Q)F  -f  0  <  0 


where  4>  = 


0  1 
1  0 


4-  = 


-1  jWc 
-juJc  -UI1UJ2 


Wc 


(a;i-|-a;2) 

2 


(9) 

(10) 


Lemma  5  [13]  Given  a  symmetric  matrix  E  G  and  two  matrices  X,  Z  of  column  dimension  p,  there 
exists  a  matrix  Y  such  that  the  LMI 

Y-y  symX'^YZ  <0  (11) 

holds  if  and  only  if  the  following  two  projection  inequalities  with  respect  to  Y  are  satisfied: 

X-^'^YX-^  <  0,  Z-^^YZ-^  <  0.  (12) 

where  X^  and  Z^  are  arbitrary  matrices  whose  columns  form  a  basis  of  the  null  spaces  of  X  and  Z,  respec¬ 
tively. 
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3.  STABILITY  ANALYSIS 

we  are  now  in  a  position  to  present  a  new  condition  for  checking  the  stability  of  2-D  continuous 
systems  of  Roesser  model. 

Lemma  6  The  2-D  continuous  system  (1)  is  asymptotically  stable  if  there  exist  Pi  >  0  and  P2  >  0  such 
that  the  LMI 

X^P  +  PA<0 


is  feasible.  Where  P  =  dtapjPi,  P2}  and  A  = 
Proof  LMI  (13)  can  be  rewritten  as 


Ai  A2 
A3  A4 


PiAi  +  Af  Pi  AfP2  +  P1A2 
*  P2A4  +  AfP2 


and  this  latter  LMI  can  be  rewritten  as 


Ai  A2 

T 

■  0  Pi  ■ 

Ai  A2 

I  0 

.Pi  0  . 

I  0 

<  0 


-f  0  <  0 


where 

0  = 

by  Lemma  3,  (15)  is  equivalent  to 


A3  A4 

T 

■  0  P2  ■ 

[As  A4  1 

0  I 

P2  0 

' — 1 

0 

_ 1 

(jwI-Ai)  1^2 

I 


0 


{jwl  -  Ai)  ^2 

I 


<0, 


(13) 


(14) 


(15) 


(16) 


'S'(jw) 

* 

0 

P2  ■ 

/ 

P2 

0 

I 

<  0 


(17) 


where  S{jw)  is  the  frequency  response  matrix  obtained  from  S{s)  of  Lemma  4,  moreover,  (17)  can  be  written 


S{jw)*P2+P2Sijw)  <0. 


(18) 


and  the  existence  of  a  P2  >0  satisfying  this  last  condition  immediately  implies  that  all  eigenvalues  of  S{juj) 
must  have  strictly  negative  real  parts,  Vw  €  M  U  00,  that  is,  feasibility  of  (13)  guarantees  that  condition  of 
Lemma  2  holds.  Moreover,  feasibility  of  (13)  implies  that 

PiAi  +  Aj^Pi  <  0,  P2A4  +  A^P2  <  0, 

and,  since  Pi  >  0  and  P2  >  0,  all  eigenvalues  of  the  matrices  Ai  and  A4  must  have  strictly  negative  real  parts, 
and  feasibility  of  (13)  guarantees  that  Lemma  1  and  Lemma  2  are  satisfied.  □ 


Lemma  6  proposes  an  LMI  condition  for  the  asymptotical  stability  of  the  system  in  (1),  there  exists  some 
conservativeness  due  to  the  requirement  of  a  constant  matrix  P2  for  all  a;  G  M  U  00,  though.  Lollowing  the 
similar  line  of  [18,  19,  20],  the  existence  of  P2{juj)  >  0  such  that 

S{jw)* P2{jw)  +  P2{jw)S{jw)  <0,  Vw  G  M  U  00, 

is  a  sufficient  condition  for  asymptotical  stability  of  the  system  (1).  Based  on  this  result,  and  in  order  to  reduce 
the  conservativeness  of  Lemma  6,  we  attempt  to  obtain  a  piecewise  constant  matrices  P2{juj)  via  a  frequency- 
partitioning  appoach,  over  the  entire  frequency  field. 

Denote  H  =  M  U  00,  and  due  to  S{—juj)  =  S{juj)*,  the  following  identities  hold: 

sup  ReX{S{jw))  =  sup  ReX{S{jw))  =  sup  ReX{S{jw))  (19) 
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where  =  [0,  oo],  =  [oo,  0].  Therefore,  it  suffices  to  consider  the  half  frequency  field  12+.  Now,  for  a 
given  positive  integer  N,  dividing  the  frequency  domain  12+  into  N  intervals,  such  that 

N 

=  =  ^0  =  0,  a;7v  =  oo,  (20) 

then  applying  the  result  of  Lemma  4  on  each  interval,  we  obtain  the  following  theorem. 


Theorem  1  For  a  given  positive  integer  N,  define  frequency  intervals  12+  as  in  (20).  System  (1)  is  asymp 
totically  stable  if  there  exist  P^j  >  0,j=l,2.  Pu,  P2i,  Qj  >  0,  2  =  1,  2, ...,  N 

A^P,  +  P,A  +  0  Qi)F  <  0 


^jxjPsj  +  PsjAjxj  <0,  j  —  1,2. 


where  A  = 


Ai  A2 
A3  A4 


F  = 


Ai  A2 

I  0 


For  i  =  2,3, N  —  1, 

For  i  =  1  and  i  =  N, 

«'i  = 


-1 

jWc 


JWcz 


and  Pi  =  diag{Pii,P2i\. 

(Wi-i+Wi) 

with  Wri  =  - z - . 


-1  0 

0  wl 


and  4'jv  = 


1 


0 


0  -w^_i 


(21) 

(22) 


(23) 


(24) 


respectively. 

Proof  By  Assumption  2,  we  have  Re\{Ai)  <  0  and  ReX{A4)  <  0  if  and  only  if  there  exist  Pgi  >  0 
and  Psi  >  0  such  that  Af  Pgi  +  PsiAi  <  0,  Aj Ps2  +  Ps2Ai  <  0,  LMIs  in  (22)  are  satisfied. 

For  i  =  2, ...,  N  —  1,  and  {i  =  1,  i  =  N},  according  to  [14]  the  matrix  T'i  should  be  taking  as  (23)  and  (24), 
respectively.  The  condition  in  (21),  can  be  written  as 

F^{^<S)Pu  +  '^>^<S)Q^)F  +  &,<0  (25) 


where 


and 


$  = 


0  1 
1  0 


0,:  = 


A3  A4 
0  / 


iT 


($  0  R 


2i) 


A3  A4 
0  I 


(26) 


(27) 


Denote  G{juj)  =  (jujl  —  Ai)  ^A2,  and  S{juj)  =  A4  +  A3(jujl  —  Ai)  ^A2  =  A3G{juj)  +  A4,  by  lemma  4, 
the  following  inequality  follows: 


■  GiM  - 

* 

0, 

■  Gijw)  - 

I 

t 

I 

5'(iw) 

I 


S(.juj) 

I 


<  0,  Vw  G  12+ 

<0,  Vw  G  12+ 


or  in  a  more  compact  form 


S{jw)*P2r  +  P2rS{juj)  <  0,  Pa*  >0,  Vw  G  12+ 


(28) 


(29) 


(30) 


So  ReX{S{ju;))  <  0  is  finally  guaranteed  for  all  w  G  12+.  Combining  PeA(Ai)  <  0,  PeA(A4)  <  0  and 
ReX{S{juj))  <  0,  we  conclude  that  system  (1)  is  asymptotically  stable  based  on  Lemma  2.  The  proof  is 
completed.  □ 
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Remark  2  When  N  =  1,  and  if  letting  Qi  =  0,  Pu  >  0  and  P2i  >  0  be  real,  then  (21)  reduces  to  (13), 
that  is.  Lemma  6  is  a  special  case  of  Theorem  1. 

Theorem  2  For  a  given  positive  integer  N,  define  frequency  intervals  as  in  (20).  System  (1)  is  asymp¬ 
totically  stable  if  there  exist  Pgj  >  0,  j=l,2.  Pu,  P2i,  Qi  >  0,  Wu,  W2i,  z  =  2, 3, -/V  —  1,  such  that 

Alii  Ai2i  0  Ai4i 

*  A22i  A23i  A24i 

*  *  Aasi  A34i 

*  *  *  A44i 

A.,h  A,,-2*1^q^  i  =  l,2.  (32) 

*  ^sjZi 

All.  =  -Q^  -  Wu  -  Wl, 

Al2i  =  Pu  +  jWciQi  —  Wu  +  WuAi, 

Al4i  =  WuA2, 

A22i  =  —Wi-lWiQi  +  A^Wu  +  W^liAi, 

A23j  =  A^Wl^, 

A24i  =  -|-  WuA.2, 

A33.  =  -VP2.  -  Wlu 

A34i  =  P2i  —  W2i  +  W2iA4, 

A44i  =  112.^4  -I-  A^W2i- 

Kjii  =  -w,i  -  wju 

Asj2i  —  Psj  ^^li  P  WjiAjxj, 

Asjsi  =  ^jxjWji  +  WjiAjxj,  j=f2. 

For  i=l,  we  replace  An^,  and  Ai2i  in  (31)  by 
A121  =  Pu  —  Wii  -I-  WiiAi, 

A221  =  wfQi  -I-  AfWu  -f-  WiiAi,  respectively. 


For  i=N,  we  replace  An^,  Ai2i  and  A22i  in  (31)  by 

Aiiiv  =  Qn  —  WiM  —  Wipf, 

Ai2N  =  -PlAT  —  WIj^  -\-  WinAu 

A22N  =  —w%_iQn  +  -I-  WinAi,  respectively. 

Proof  From  Theorem  1,  let 


$  0  Pi.  +  4-.  0  Qi  0 

*  $  0  P2i 


According  to  [14],  for  z  =  2, A^  —  1,  i  =  1  and  i  =  N,'^i  as  in  (23),  (24),  and  $  as  in  (26). 
■  11^1.  0  1 

T  ll^li  0  ry  — Al  0  A2  T  /->  1  \  •  "I 

Let  F  =  „  ...  ,Z=  „  .  T  A  ,  X  =  I,  (31)  IS  equivalent  to 

0  11^2.  [  0  A3  -/  A4  J’  ^ 

0  1F2. 


,  X  =  I,  (31)  is  equivalent  to 


sym{X^YZ)  +  S  <  0  (34) 

since  one  can  choose  X-^  =  0,  the  first  inequality  in  (12)  vanishes,  and  then  by  lemma  5,  (34)  hold  for  some  Y 


if  and  only  if  Z^  YZ^  <  0.  Note  that  Z^  can  be  selected  as  = 


Al  A2 
1  0 
A3  A4 
0  7 


,  and  then  by  calculation,  we 


can  obtain  the  equivalence  between  YZ^  <  0  and  (21).  Consequently  (21)  is  equivalent  to  (31). 
In  addition  from  (22),  we  get 


0  P 


'a  3  1  r  Ajxj 


Paj  0 


<0,  j  =  l,2. 
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Figure  1:  ReXmax{S{juj))  and  7* 


The  equivalence  between  (35)  and  (32) 

E  = 

0  Psj 

,  Y  = 

Psj  0 

Wji  _ 

be  similarly  found  by  re-introducing 
Z  =[-I  ],  X  =  I,  j=l,2. 


Theorem  1. 


Thus,  Theorem  2  is  equivalent  to 

□ 


Example  1  In  this  part,  we  provide  an  example  to  illustrate  the  application  of  the  proposed  method, 
consider  system  in  (1),  where  the  matrices  in  the  system  are  obtained  by  a  suitable  transformation  from  the  original  system 
matrices  [21  ]: 

t  (-4- 

The  matrices  in  the  original  problem  are  as  the  following  form  [18]: 


A,  = 


Adi 

Ad2 

■  0.5 

0.5 

■  0.4 

1.1  ■ 

■  -0.1 

-0.1  ■ 

Ad  = 

Ads 

AdA 

,  Adi  ~ 

0.1 

-0.1 

,  Ad2  = 

0.6 

0.1 

,  Ads  = 

-0.2 

0.6 

Adi  = 


-0.5  -0.5 

-0.1  -0.7 


we  obtain 


= 


Ai  — 


5.2979  -16.0426 

4.2128  -12.7447 

-22.5745  -23.4894 

-26.9787  -30.5745 


A2  = 


20.5957 

16.4255 


23.9149 

16.5106 


As  = 


-5.7872 

-7.4894 


16.2553 

22.2128 


Denote  7*  as  the  minimum  value  of  —ji  that  satisfies 


sup  ReX„ 


<  -7i  <  0 


7*  could  be  computed  from  (31)  by  replacing  >1>  (g)  Psi  in  (33)  by 


0  P2. 

Psi  Ji 


and  minimizing  —Ji-  Figure  1  shows 


ReXmax{S{ju}))  and  the  executed  7*  by  Theorem  2  with  N  =  1,2, 4,  8.  The  stability  of  the  above  system  is  verified,  since 
Re{Xmax{S{ju:)))  <  0  is  evident.  With  N  growing,  It  is  further  shown  that  —7;  tends  to  the  value  of  ReXmax{S{juj)) 
over  TT^. 


Theorem  2  with  N  =  1  fails  to  decide  the  stability  of  the  above  system.  By  increasing  N ,  it  is  found  that  Theorem 
2  with  N  =  2, 4,  8  succeeds,  note  that,  the  above  system  is  asymptotically  stable  only  for  i  =  2,  ...,N.  But  for  i  =  1, 
whatever  N,  and  whatever  the  way  of  partitioning  the  entire  interval,  always  system  is  not  stable.  This  is  due  to  the  rapid 
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variation  of  the  curve  of  Re\max{S{jiS))  in  the  vicinity  ofujo  =  0. 


Remark  3  in  the  following  domain  [16,+oo],  ReXmax{S{juj))  remains  relatively  stationary  to  the  value  ReXmax{S{jc 
ReXmaxiA4)  =  —1.0850,  then  7*  also  tends  to  this  value  throughout  the  domain.  Even  if  it  decomposed,  we  find  very 
similar  values  to  —1.0850.  That’s  why  we  worked  on  just  the  domain  [0, 16]  (Figure  1 ). 


CONTROL  LAW  DESIGN 

In  this  section.  Theorem  2  is  further  developed  for  state-feedback  control  of  the  uncertain  2-D  continuous  sys¬ 


tems. 

Consider  a  2-D  continuous  system  of  Roesser  model  with  norm-bounded  uncertainty: 

-b 


dx''{ti,t2) 

dti 

dt2 


A.I  +  A2  ^A.2 

A3  +  AA3  A4  +  AA4 


a: '“(11,12) 
*”(11,12) 


Bi  +  ABi 

B2  A  AB2 


m(1i,12) 


where  the  uncertain  matrices  AA,,  g  =  1,  2,  3, 4  and  ABp,  p  —  1,2  formed  as 


[  AAi  AA2  ABi  ]  =  HiA  [  El  E2  Li  ] 
[  AA3  AA4  AB2  ]  =  H2A  [  E3  E4  L2  ] 


(36) 


(37) 


where  Hi,  H2,  Ei,  E2,  E3,  E4,  Li  and  L2  are  known  constant  matrices,  A  is  norm-bounded  parameter  uncertainty  satis¬ 
fying  A^  A  <  I.  Suppose  the  system  (36)  is  controlled  by  a  state-feedback  controller: 


1(11,12)  =[Ai  K2] 


*'“(11,12) 

*”(ll,l2) 


where  Ki  and  K2  are  the  controller  gains  to  be  found,  then  the  closed-loop  system  is  given  by: 


dx'^(ti,t2) 

dti 

dx"^  (ti,t2) 
dt2 


A.cl  “h  Ac2  ^Ac2 

Ac3  -h  A^c3  Ac4  +  AAc4 


*'“(11,12) 

*”(ll,l2) 


(38) 


(39) 


where 

Aci  =  Ai  A  BiKi,  Ac2  =  A2  A  B1K2,  Ac3  =  A3  -I-  B2K1,  Ac4  =  A4  -I-  B2K2, 

AAci  =  AAi  A  ABiEi,  AAc2  =  AA2  -t-  AB1E2,  AAcs  =  AA3  -t-  AB2E1,  AAc4  =  AA4  -t-  AB2H2. 

Our  objective  is  to  find  a  state-feedback  controller  in  (38)  for  the  system  (36)  such  that  the  closed-loop  system  (39)  is 
asymptotically  stable  for  all  possible  uncertainties.  Before  we  proceed,  the  following  lemma  which  is  usually  used  in  the 
robust  control  of  systems  will  be  given  first. 

Lemmal  [22]  Let  T,i,T,2  and  Abe  real  matrices  with  appropriate  dimensions  such  that  A^  A  <  I.  Then,forany 
scalar  e  >  0  the  following  inequality  holds: 

S1AS2  A  Ei2  a  Sx  <  s  XjiXji  -|-  s^2  ^2  (40) 

Now,  based  on  Theorem  2,  we  have  the  following  analysis  result  on  robust  stabilization  of  the  2-D  continuous  system  (39). 

Proposition  1  For  a  given  positive  integer  N,  define  frequency  intervals  fl'i  as  in  (20).  System  (39)  is  asymptoti¬ 
cally  stable  for  all  A  satisfying  A'^A  <  I,  if  there  exist  matrices  Psj  >  0,  y  =  1,  2.  Ps2  >  0,  Pu,  P2i,  Qi  >  0,  Wi,  W2, 
and  scalars  ei  >  0,  i  =  1,2, ...,  N,  such  that 


—1  —2 

*  H3 


— sjl  “Sj2 

*  Ssj3 


<  0 

<0,  i  =  l,2. 


(41) 

(42) 


All 

Ai2 

0 

Ai4 

* 

A22 

A23 

A24 

_ 

* 

* 

A33 

A34 

1  “2 

* 

* 

* 

A44 

Qi  - 

Wi  - 

ITf. 

li  +  j^ciQi 

-  Wf 

-b  WiA,i, 

H2  =  [  Ti  eiT 


iT2  1,  S3  =  diag{—£il, —Sil, —Eil , —£il}. 


Ai4  =  IT'iAc2, 

A22  =  —oJi-iUiQi  A  A^iW^  -b  ILiAci, 


Stability  and  Robust  Stabilization  of  2-D  Roesser  Continuous  Systems  ...  (Ismail  Er  Rachid) 


998 


□ 


ISSN:  2088-8694 


’  WiA,2, 


A23  =  A^3W2^ 

A24  =  A^3W2^ 

A33  =  -W2-  W2^, 

A34  =  P2i  —  +  11^2  Ac4, 

A44  =  W2A,4  +  A'^iWi'. 


■  WiHi 

0 

0 

0 

Ti  = 

WiHi 

0 

0 

W2H2 

,  T2  = 

(EjPLjK^)'^ 

0 

(E3PL2Kj)'^ 

0 

0 

W2H2 

(E2  -h  LiA'2)^ 

(Ej  P  L2A'2)^ 

“Sjl 

Asji 

Asj2 

Asj3 


A, 


sjl 

* 


sj2 


sj3 

T 


H.,2  =  [  T, 


sjl  sj2  ])  “sj3  —  (iicig{  Sil  ^  £i/}. 


=  -W2  -  W,  , 

=  Psj  -Wl  + 

=  Al^jWj  +  W,A, 


Ts21  = 


W,H, 


Ts22  = 


(Ejxj  +  LjKj) 


=  1,2. 


Proof  From  Theorem  2,  by  replacing  VFii  and  W2i  with  FFi  and  IF2  respectively.  System  (39)  is  asymptotically 
stable  if  there  exist  Psj  >  0,  j  =  1,2.  Pu,  P2i,  Qi  >  0  such  that  inequalities  in  (3 1)  and  (32)  satisfied,  in  which  Aq  should 


be  Acq  +  AAcq  for  q  =  1,  2,  3, 4.  The  above  LMIs  can  be  re-written  into  the  following  form: 

Hi  +  TiATl’  -f  T2A^Tf  <  0.  (43) 

H,,i  -h  T«,-i ATJ,.2  -f  T,^2 A^tJ,-!  <0.  j  =  1,  2.  (44) 

According  to  Lemma  7,  the  above  inequalities  holds  for  all  A  if  and  only  if  there  exist  some  scalars  ei  >  0  such  that 

Hi  -h  e-^TiTr  -h  eiT2T^  <  0  (45) 

“sjl  P  £i  PsjlPsjl  p  £i'^ sj2P' Sj2  <  0  (46) 

which,  by  the  Schur  complement  in  [23],  (45)  and  (46)  give  rise  to  (41)  and  (42).  □ 


Retnurk  4  it  is  interesting  to  note  that,  as  the  LMIs  including  their  proofs  for  i  =  1  and  i  =  N  are  similar  to  those 
forcasesofi  =  2,...,N  —  1,  we  give  these  Proposition  I  in  one  unified  form  for  all  possible  value  of  if  or  reason  of  space. 
The  same  expression  applies  to  the  following  Theorem  3. 

Now,  based  on  Proposition  1 ,  we  are  in  a  position  to  give  a  new  method  of  state-feedback  stabilization  controller  design  for 
the  Reosser  model. 

Theorem  3  For  a  given  positive  integer  N,  define  frequency  intervals  fl’*'  as  in  (20).  System  (39)  is  asymptotically 
stable  for  all  A  satisfying  A^  A  <  I,  by  a  state  feedback  controller  in  (38),  if  there  exist  matrices  Psj  >  0,  j  =  1,  2.  Pu, 
P2i,  Qi  >  0,  Wi,  W2,  Ni,  N2  and  scalars  Si  >  0,  i  =  1,  2, ...,  N,  such  that 


(47) 


“sj 


^sj2 

“sj3 


<  0 


All 

Ai2 

0 

Ai4 

* 

A22 

A23 

A24 

“1 

* 

* 

A33 

A34 

* 

* 

* 

A44 

All  =  - 

-Qi  - 

Wl  - 

Wi^, 

Ai2  =  Pli  P  joJcQi  —  Wl  -|-  AlWj  +  BjNl, 

Ai4  =  A2Wf  p  B1N2, 

A22  =  — oJi—iiOjQi  P  IViAj^  P  AilP^  P  i?iAli  -f  JV^ PT ’ 
A23  =  P  Nf  Bj, 

A24  =  A2Wf  PWiAj  P  NfB^  p  B1N2, 

A33  =  -W2  yJVi',  _ 

A34  =  P2i  —  IV2  -|“  AjlV'^  -f  B2N2, 


H3  =  diag{—5il ,  —Sil,  —Sil,  —Sil} 


(48) 
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Figure  2:  Closed-loop  responses  of  and  Xi{ti,t2)- 


A44  =  A4Wi'  +  B2N2  -f  W2AI  +  N^B2. 

0 


Ti  = 


'  Hi 

0 

Hi 

0 

0 

H2 

.  T2  = 

0 

H2 

0 


Wi  Ef  -K  Nl  Lf  Wi  Ej  +  NIlT 


^  0 
WiE^  +  N^Ll 


0 


^sjl 


A, 


sjl 

* 


H.42  =1 


[  SE 


sjl 


T 


Asj2 

Aaj3 

Asji  =  -IF4 
Asj2  ~  Psj  ~  Wj  -f  Aj 

A., 3  =  W,AJ^^  +  Aj^.Wf  +  BjN,  +  NfBj 


sj2 


W2EJ;  +  L'^}  J 

j,  Esj3  =  diag{—6il,  —Sil}, 


^sj3 


,IF/  +BjNj, 

■J  , 


r  H, 

0 

J 

.  Hi  . 

,  Psi2  = 

_  WjEj^j  +  Nf Lj  _ 

i  =  l,2. 


If  the  above  conditions  are  satisfied,  a  stabilizing  control  law  [  K\  K2  ]  is  given  by 
Ki  =  NiWf'^,  K2  =  N2Wf'^. 


Proof  If  (42)  holds,  IFi  and  W2  are  nonsingular.  Pre-and  post-multiplying  (41)  by  nonsingular  matrices: 

diag{Wr\Wr\  Wi\  anddiag{Wr^,  Wf^,  Wf^, 


and  Pre-and  post-multiplying  (42)  by  nonsingular  matrices: 

diag{W-\W-\s-^I,s-^I}  and  diag{W-^,W-^,e-^I,e-^I},  j  =  1,2. 

Making  change  of  variables  as  follows: 

Wi  =  =  Si  =  ef\  Q,  =  WiQ.WT,  Pu  =  WiPuWP,  P2^  =  W2P2iWf , 

Psl  =  WlPs^WT,  Ps2  =  W2Ps2Wf 


We  can  obtain  the  equivalence  between  Theorem  3  and  Proposition  1,  where 

Ni  =  KfWp,  N2  =  K2Wf.  □ 


In  the  following,  we  provide  an  example  to  demonstrate  the  effectiveness  of  the  proposed  method  in  this  section. 


Example  2  Consider  the  uncertain  2-D  continuous  system  in  (39)  with  parameters  given  by:  [11] 


■  -1.2 

0.3 

-0.7  ■ 

■  -0.7 

0.2  ■ 

= 

-1 

0.5 

0.6 

,  A2  = 

0.5 

1 

0 

0.2 

-1.8  _ 

_  -0.5 

0 

'  -0.3 

-0.1 

0.5  ' 

1 

0.3 

-0.6 

Bi  = 

1 

1 

0.5 

0.2 

0 

0.6 

,B2  = 

—  1 

1 

El  = 

0.1  0 

0.2 

],£;2  = 

=  [0.1 

-0.2  ] 

Li  = 

,  A3  = 


0.2 

0.5 


0.9 

0 


0 

0.2 


-1.5 

0.1 


,A4  = 


-0.8 

0.1 


0.1 

,Hi  = 

0 

,  H2  = 

■  0.1  ■ 
0 

-* 

-0.2 

0  1,  E3  =  El,  E4  =  E2,  L2  = 


0.2 

0.6 


Li. 


Because  the  eigenvalues  of  matrices  Ai  and  A4  contain  positive  eigenvalues  given  by  0.4072  and  0.6141,  re¬ 
spectively.  Therefore,  the  nominal  2D  continuous  system  under  consideration  is  not  asymptotically  stable.  The  aim  of 
this  example  is  to  design  a  frequency-partitioning  state  feed-back  controller  such  that  the  resulting  closed-loop  system  is 
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asymptotically  stable  for  all  admissible  uncertainties.  By  application  of  Theorem  3  on  2-D  contunuous  Roesser  model,  we 
obtain  the  state  feedback  controller  parameters  as  results: 


■  0.3643 

-0.2418 

-0.5541  ■ 

'  -0.1881 

-0.5412  ■ 

0.1982 

-0.2928 

1.0384 

,  K2  = 

0.3719 

1.3484 

0.2830 

-0.6519 

2.4311 

0.1270 

-1.8949 

To  show  the  asymptotic  stability  via  the  state-feedback  control  in  (38)  with  (47,48),  the  state  evolution  ofx^{ti,t2)  and 
x'’{ti,t2)  of  the  closed-loop  system  in  (39)  are  depicted  in  Figure  2.  For  simulations,  assume  that  12)  =  0  and  let 
the  initial  and  boundary  conditions  to  be: 

a;'*(0,t2)  =  0.2;  0  <  ta  <  10 

a;„(ti,0)  =  0.2;  0  <  ti  <  10 

x'^(0,t2)  =  Xv{ti,0)  =  O',  ii,t2  >  10. 

The  simulation  results  show  that  the  closed-loop  system  in  (39)  is  asymptotically  stable. 


5.  CONCLUSIONS 

This  paper  has  been  tackled  the  stability  and  robust  stabilization  problem  of  2-D  continuous  Roesser  systems. 
The  proposed  conditions  of  the  system’s  stability  has  been  provided  in  terms  of  LMIs.  Furthermore,  GKYP  lemma  com¬ 
bined  with  frequency-partitioning  approach  is  used  to  reduce  the  conservativeness.  Robust  stabilization  using  static  state 
feedback  has  been  studied  as  well,  and  the  stabilizing  feedback  gain  matrices  have  been  constructed  based  on  the  solutions 
of  certain  LMIs.  Finally,  numerical  examples  demonstrate  that  the  proposed  methods  are  effective. 
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